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ABSTRACT 

Fast particles diffusing along magnetic field lines in a turbulent plasma can diffuse 
through and then return to the same eddy many times before the eddy is randomized 
in the turbulent flow. This leads to an enhancement of particle acceleration by large- 
scale compressible turbulence relative to previous estimates in which isotropic particle 
diffusion is assumed. 



1. Introduction 

Particle acceleration by magnetohydrodynamic (MHD) turbulence has been investigated by a 
number of authors (e.g., Fermi 1949, Hall & Sturrock 1967, Wentzel 1968, Kulsrud & Ferrari 1971, 
Eichler 1979, Achterberg 1981, Schlickeiser & Miller 1998), and has been studied in a number 
of astrophysical settings, including the interstellar medium (Ptuskin 1988, Schlickeiser & Miller 
1998), solar flares (LaRosa et al 1994, Miller et al 1996, Miller et al 1997, Miller 1998, Schlickeiser 
& Miller 1998, Miller 1998), active galactic nuclei (Blackman 1999, Gruzinov & Quataert 1999), 
diffuse intracluster plasma in galaxy clusters (Eilek & Weatherall 1999, Weatherall & Eilek 1999, 
Dogiel 1999, Brunetti et al 1999, Ohno et al 2002, Fujita et al 2003), and the lobes of extragalactic 
radio sources (Eilek & Henriksen 1984, Manolakou et al 1999). 

The way in which magnetohydrodynamic (MHD) turbulence accelerates energetic particles 
depends on the relative magnitudes of and X^fp, where is the length of a turbulent eddy (or 
the wave length of a weakly damped wave) measured along the magnetic field, and Xmfp is the 
energetic-particle scattering mean free path. In the weak-scattering limit, ^ ^imfp^ a particle's 
velocity V|[ along the background magnetic field is approximately constant as the particle travels 
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many wave lengths along the magnetic field (assuming that the amplitudes of the fluctuations are 
sufficiently small), and the particle interacts resonantly with a weakly-damped wave when the 
wave frequency in the frame of reference moving along the magnetic field at speed V|| is an integer 
multiple of the particle's gyrofrequency. In the strong-scattering limit, ^ Knfp-> ^ particle travels 
diffusively across an eddy of length Xy, and acceleration is non-resonant. 

Particle acceleration in the strong-scattering limit is the focus of this paper. Ptuskin (1988) 
studied particle acceleration in the strong- scattering limit by a power-law spectrum of acoustic 
or magnetosonic waves extending from a large scale / to a smaller scale li , with rms fluctuation 
velocity Wrms dominated by the waves at scale /. He assumed that particle diffusion in space is 
isotropic with diffusion coefficient D, as opposed to primarily along field lines, and that some 
type of small-scale wave distinct from the acoustic/magnetosonic waves scatters the particles, with 
^mfp ~ Z)/v <S /i, where v is the particle speed. He found that for D » Vw/, where is the wave 
phase velocity, particles diffuse in momentum space according to the equation 

^ = lAf/z) ^ 

dt dp \ ''dp 



P^D' ] , (1) 



where / is the particle distribution function averaged over the waves, and the momentum diffusion 
coefficient Dp is given by 

n^^P^ (2) 

where p is the particle momentum. The condition D^Vy,l implies that particles diffuse a distance / 
in a time much shorter than the wave period at scale /. In this limit, the largest-scale waves make 
the dominant contribution to Dp. 

This paper extends Ruskin's analysis to strong compressible MHD turbulence taking into 
account the reversible motion of particles along magnetic field lines. It is assumed that the rms 
turbulent velocity Wrms equals or exceeds the Alfven speed va, and that there is a single dominant 
length scale (outer scale) denoted / for both velocity and magnetic fluctuations. Only the fluctu- 
ations at scale ~ / are considered. As in Ptuskin's (1988) study, it is assumed that some type of 
small-scale wave scatters particles, and that Xmfp ^ I- It is assumed that eddies are randomized 
in a time Xrand ~ V"rms, and that Dp ^ Urmsh where D|| is the diffusion coefficient for particle 
motion along field lines. Thus, a particle diffuses a distance greater than / during a time Xrand^ and 
can diffuse through and then return to an eddy multiple times before the eddy is randomized in 
the turbulent flow. It is shown in section 2 that this increases the coherence of stochastic particle 
acceleration, effectively increasing Dp above the value in equation (2) by a factor of ~ A^, where 

~ ^D^^/ul is the typical number of times a particle diffuses through and then returns to an eddy 
before the eddy is randomized, assuming a particle is tied to a single field line. The discussion of 
section 2 consists of estimates based on physical arguments; these estimates are tested and vali- 
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dated by Monte Carlo particle simulations in section 3. In section 4 the results of the paper are 
applied to particle acceleration in the interstellar medium and the lobes of radio galaxies. 

It is known that small cross-field motions and the divergence of neighboring magnetic field 
lines in turbulent plasmas allow a particle to escape from its initial field line by traveling sufficiently 
far along the magnetic field. This places an upper limit on N, which is explained in the Appendix. It 
is shown that this upper limit is not important in the astrophysical examples considered in section 4 
given the assumed values of However, for sufficiently large values of D||, the upper limit on N 
becomes important. 



2. Phenomenological description of momentum diffusion 

Since the scattering mean free path is assumed to be <^ / and only those eddies at scale ~ / 
are considered, a particle's momentum p evolves according to the equation 

where 

y = Hp/po), (4) 

po is a particle's initial momentum, and u is the turbulent velocity associated with the large-scale 
eddies (Skilling 1975, Ptuskin 1988). As a particle diffuses in space, it encounters a series of 
eddies that cause random increments iny, leading to diffusion in y-space with diffusion coefficient 



2 

rms 



where dt is the time during which increments to y remain correlated, and (63;)rms is the rms change 
in y during a time dt. It is assumed that the turbulent velocities are randomized in a time 

bl 

^rand = , (6) 

Mrms 

where & is a constant of order unity. The characteristic time for a particle to diffuse through an 
eddy is 

It is assumed that Dm » Mrms^ and thus 



■Tdiff -C 'Crand- 



(8) 
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In this section, it is assumed that a particle remains tied to the same magnetic field line. Thus, 
given equation (8), a particle typically diffuses through an eddy and later returns to the same eddy 
before the eddy is randomized. During such successive encounters, the values of V ■ u within the 
eddy are correlated, and the corresponding increments to y are correlated. The correlation time dt 
for random steps in _y-space is thus Xrand> and not Xdiff. 

The increment to y during a time Trand is 

1 /-trand 

5y=-T/ V-uJf. (9) 
3 Jo 

It is assumed that there is no average compression or decompression in the turbulent flow, so that 

(53;) = 0, (10) 

where (...) indicates an average over an arbitrarily large number of random steps. During a 
time Trand. ^ particle's position is distributed over an interval of length ~ Nl, where 



Ar=^5^. (11) 

Since the values of V • u in different eddies are uncorrelated, the value of V • u averaged over a 
distance Nl is ~ iV-i/2(V-u)nns, where (V-u) rms is the rms value of V • u. Thus, the rms value of 
8y is approximately 

A^"^/^(V-u)nnsi;rand ,1^, 
(Oj)rms ^ • (12) 

Since only eddies of scale / are considered, and since Mrms ^ va, 

(V-U)nns=— (13) 

where a is a constant of order unity. Thus, ifN:^ 1, then (5j)rnis ^ 1, and the probability P{y, t)dy 
that ln(p/ pq) e {y-iy + dy) at time t satisfies the Fokker-Planck equation 

dP 3^ 

where Dy, which is defined in equation (5), is given by 

Dy^"^^. (15) 

^ 18Z)|| ^ ^ 

The log-momentum diffusion coefficient Dy is analogous to Ptuskin's (1988) momentum dif- 
fusion coefficient £)p,Ruskin divided by p^. If Ruskin's coefficient D for isotropic spatial diffusion 
is set equal to Dm, then in the large-A^ limit 



iV. (16) 

^p,Ruskin 



P^Dy 
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This ratio can be understood as follows. Ruskin's result is approximately recovered by taking 
dt ~ Xdiff. Such a relation is valid if spatial diffusion is isotropic, since there is only a small chance 
that a particle returns to the same eddy multiple times before the eddy is randomized when D ^ 
/t/rms- However, when a particle diffuses along a single field line, the particle's location during a 
time Xrand = ^''^'^diff is distributed among roughly eddies, and thus a particle spends a time ~ 
NXdiff in an eddy before the eddy is randomized. The particle's random walk in y-space can thus be 
thought of consisting of individual steps of duration NXtUff and magnitude dy — A'^Xdiff (V ■ u)rms/3. 
The resulting value of Dy is a factor of larger than in the isotropic-spatial-diffusion case in which 
a particle takes individual random steps in y-space of duration Xdiff and magnitude 8y = Xdiff (V • 
u)nns/3. Thus, the reversible motion of a particle along a field line makes the particle's random 
walk in y-space more coherent, and its diffusion in y-space more rapid, than in the isotropic-spatial- 
diffusion case. 

One consequence of equation (15) is that — > oo as Xrand °°- The reason is that if V • u were 
constant in time along a field line, increments to y would remain correlated forever if a particle were 
tied to that field line, and the particle would undergo superdiffusion in j- space. This phenomenon 
can be described approximately as follows. During a time t, a particle's position is distributed 
over ~ M eddies, where 

M=-^^j^. (17) 

Since the values of V • u within different eddies are uncorrected, the typical value of V • u averaged 
over M eddies is ~ M~^/^(V • u)nns- For large M, the rms value of y for a particle obeys the 
equation 

dyrms M-l/2(V-u)nns ,-1/4 

dt 3 
Integrating and squaring the result, one obtains 



o^t-'i\ (18) 



,2 1 6/(V-u)L?3/2 



3. Numerical simulations 

In this section, particle acceleration by compressible turbulence is investigated with the use 
of Monte Carlo particle simulations. Particles are assumed to follow a single field line. Distance 
along the field line is denoted by the coordinate x. The values of V • u at space-time grid points 
(x,t) = (m/,nXrand) ^rc set to either +x~^j or — x"^^ with equal probability, where m and n are 
integers. The value of V • u between grid points is obtained by linear interpolation in space and 
time, which gives (V • u)nns'Crand = 2/3, or, equivalently, ab = 2/3, where a and b are defined in 
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equations (6) and (13). During each time step of duration bt, a pari;icle's y- value is incremented by 
an amount —V • u5?/3, and there is an equal probability that its x-coordinate is either increased or 



decreased by an amount 




In figure 1, Z)y is plotted for different values of Dy. To express the results in terms of Wrms 
and /, it is assumed that a = 1. The values of Dy in the numerical simulations, depicted with 
solid triangles in figure 1, are obtained by dividing the average of for 10"^ simulated particles 
by 2t att = IBOXrand- The solid line represents the value of Dy in equation (15). The dashed line 
represents the scaling Dy analogous to Ptuskin's (1988) scaling Dp oc D~^ for isotropic 

spatial diffusion of particles. 




Fig. 1 . — Dy as a function of . The solid triangles are the results of the Monte Carlo simulations. 
The solid line represents equation (15), and the dashed line represents a scaling of Dy oc D'j^^. 

The increase of y^^^ with time when V ■ u is constant in time along a field line is shown in 
figure 2. The solid line is the value of y^^^ in equation (19). The averages of y^ for 10^ simu- 
lated particles are given by the solid triangles and show superdiffusive behavior in approximate 
agreement with equation (19). 
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Fig. 2. — Superdiffusion in log-momentum-space when V ■ u is static along a field line. The solid 
triangles give the average of for 10^ simulated particles that are tied to a single magnetic field 
line. The solid line is a plot of equation (19), in which oc ?3/2 



4. Application to Galactic cosmic rays and relativistic electrons in radio galaxies 

From equations (14) and (15) the characteristic time for particle acceleration by large-scale 
compressible turbulence is 

18d!/2/V2 

Tacc , (20) 

Mrms 

where it is assumed that a = 1 in equation (13) and b = I in equation (6). The values of Z)|| for 
Galactic cosmic rays and relativistic electrons in the lobes of radio galaxies are not well known. 
There is disagreement as to the correct way to model the waves that cause energetic particle scat- 
tering (see, e.g., Bieber et al 1994, Schlickeiser & Miller 1998, Chandran 2000, Felice & Kulsrud 
2001, Yan & Lazarian 2002). Some studies indicate that D and Dp have different values in differ- 
ent phases of the ISM (e.g., Felice & Kulsrud 2001), and the values of D and D|| may be different 
in the disk and halo of the Galaxy (Berezinskii et al 1990). In the hot coronal plasma of the 
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ISM, Mrms ~ 10^ cm/s (Shelton et al 2001). This value is comparable to the sound speed, and 
thus the turbulence is compressible. It is assumed that / ~ 3 x 10^*^ cm. If D|| = 10^^ cm^/s, then 
'Cacc = 3 X 10^ yr. If Z)|| = 10^^ cm^/ s, then Xacc = 10^ yr. These values of Xacc are comparable to the 
Galactic confinement time of 1-GeV cosmic rays, and are the types of time scales required if reac- 
celeration is to explain the peaking of secondary/primary ratios at a few GeV (Ptuskin 2001). For 
relativistic electrons in radio lobes, it is simply assumed that = VA/iobe^ where kobe — 100 kpc 
is the typical size of a radio lobe. For this value of D||, the time for electrons to diffuse out of a 
radio lobe is of order the Alfven crossing time of the lobe. The thermal plasma density is taken to 
be 5 X 10~^ cm~^, which is an an upper limit based upon Faraday depolarization, and B is taken to 
be 20 //G, a typical equipartition field (Spangler & Pogge 1984, Spangler & Sakurai 1985). These 
values yield va = 6 x 10^ cm/s, and it is assumed that Mrms = va- The value of / for radio-lobe 
turbulence is taken to be /iobe/30 based on the standard deviation of the linear Stokes parameters in 
3C 166 (Spangler 1983). For these parameters, Xacc ~ 5 x 10^ yr. This type of time scale is com- 
parable to the ages of relativistic-electron populations obtained from synchrotron aging estimates, 
which suggests that reacceleration may have to be taken into account in models of radio-lobe elec- 
trons. 



5. Summary 

It is shown that the reversible motion of energetic particles along the magnetic field enhances 
the rate of particle acceleration by large-scale compressible turbulence relative to the isotropic- 
spatial-diffusion case. When this effect is taken into account, the time scale for acceleration 
of relativistic electrons by large-scale compressible turbulence in the lobes of radio galaxies is 
~ 5 x 10^ yr, and the time scale for acceleration of 1-GeV protons by large-scale compressible 
turbulence in the coronal plasma of the ISM is ~ 3 - 10 X 10^ yr. The main source of uncertainty 
in these estimates is the uncertainty in ; other sources include the uncertainty in /, Wrms, and the 
dimensionless constants a and b that relate Xrand and (V • u)nns to Mrms and /. 

I thank Eliot Quataert, Steve Spangler, Torsten Ensslin, and Eric Blackman for helpful com- 
ments. This work was supported by NSF grant AST-0098086 and DOE grants DE-FG02-01ER54658 
and DE-FC02-01ER54651. 



A. Particle escape from field lines through cross-field motion and field-line divergence 

Particles with gyroradii much less than the dominant magnetic field scale length / drift 
slowly across magnetic field lines because of gradients in the magnetic-field strength, curvature 
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of field lines, and wave pitch-angle scattering. Cross-field motions enhanced by the divergence of 
neighboring field lines allow particles to "escape" from their initial field lines, that is, to separate 
from their initial field lines by a distance I (Rechester & Rosenbluth 1978, Chuvilgin & Ptuskin 
1996, Chandran & Cowley 1998). This can be seen with the aid of figure 3. Suppose a particle 




Fig. 3. — Diagram showing the way in which the divergence of neighboring field lines allows a 
particle to escape from its initial field line. 

starts out at point P on field line Fi, and that its diffusive motion along the field initially moves 
it towards point Q. After moving a short distance, field gradients and collisions cause the particle 
to take a step of length ~ Pg across the magnetic field, from Fi to a new field line F2. Although 
the particle continuously drifts across the field to new field lines, let us assume for the moment 
that it remains attached to F2. As the particle diffuses along F2, F2 diverges from Fi. Let Zs be 
the distance that F2 must be followed before F2 separates from Fi by a distance /. (Because the 
particle continuously drifts across the field, it typically separates from Fi after traveling a distance 
somewhat less than Zs along the field; this effect is ignored in this discussion.) After the particle 
moves a distance Zs along F2, its subsequent motion is not correlated with Fi. The particle proceeds 
to point R, and then its diffusive motion along the magnetic field changes direction, bringing it 
back towards point Q. Instead of following F2 back to point Q, however, the particle drifts across 
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the field and ends up on a new field line F3. After following F3 for a distance ~ Zs, the particle 
separates from F2 by a distance ~ / and proceeds to point S. In this example, after the particle 
moves significantly farther than ~ Zs along the magnetic field, it has a negligible probability of 
returning to within a distance ~ / of its initial location. The maximum time a particle can remain 
correlated with an eddy is thus of order 

' - ^^'^^ (Al) 



where (zs) is the average value of Zs over a large sample of field-line pairs. 

Several studies have calculated (zs)- Jokipii (1973) and Skilling, Mclvor, & Holmes (1974) 
found that (zs) ~ I for isotropic MHD turbulence with SB > Bq, where 55 is the fluctuating mag- 
netic field and Bq is the strength of any background field in the system that is coherent over dis- 
tances » /. Narayan & Medvedev (2001) found {zs) ~ / for locally anisotropic MHD turbulence 
with 8B Bq. a related quantity, the Kolmogorov entropy, which describes the exponential diver- 
gence of neighboring field lines while their separation is smaller than the dissipation scale of the 
turbulence, has been calculated by several authors (e.g., Rechester & Rosenbluth 1979, Zimbardo 
et al 1995, Barghouty & Jokipii 1996, Casse, Lemoine, & Pelletier 2001). Another related quan- 
tity, the diffusion coefficient for the wandering of single magnetic field line, has been studied by a 
number of authors (e.g., Matthaeus et al 1995, Michalek & Ostrowski 1998). 

Maron, Chandran, & Blackman (2003) and Chandran & Maron (2003) calculated {zs) by 
tracking field lines in direct numerical simulations of incompressible MHD turbulence with an 
inertial range extending from a dominant length scale / to a much smaller dissipation scale Ij. In 
simulations with Bq = 0, they found that (zs) asymptotes to a value of order several / as is 
decreased towards the dissipation scale in the large-/ / la limit. From figures 5 and 7 of Chandran 
& Maron (2003), this value is ~ 5 — 11. We assume that {zs) has a similar value in turbulence with 
a mean magnetic field provided 8B>Bq. 

The escape of particles from field lines modifies the estimate of Dy presented in section 2 
when iy/)||Xrand > (zs)- In this case, a particle remains correlated with a magnetic field line for a 
time Tesc that is less than Xrand- During the time Xesc the particle travels a distance ~ (zs), and its 
location is spread out over roughly 

N'=^ (A2) 

eddies. The arguments leading to equations (14) and (15) can be repeated, replacing Xrand with Xesc, 
to obtain 

^^^^sdT- ^^^^ 
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To summarize, if Tesc < T^rand. then equation (A3) is valid. If Tesc > T^rand. then equation (15) is 
valid. 

For the astrophysical examples considered in section 4, Xgsc > t^rand^ and thus equation (15) 
can be applied. However, for larger values of Dy, Tesc can be reduced below Trand> in which case 
cross-field motion becomes important and equation (A3) applies. 
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